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1.  INTRODUCTION 

In  surface  acoustic  wave  (SAW)  devices,  the  interaction  between 
laser  beams  and  acoustic  waves  must  be  understood  by  considering  the 
allowed  surface  modes  of  both  disturbances.  In  this  report,  the 
formalism  for  determining  the  modes  for  electromagnetic  (laser)  waves  is 
reviewed  and  solved  for  propagation  in  an  anisotropic  dielectric  slab. 
It  is  shown  that  pure  transverse  electric  (TE)  and  transverse  magnetic 
(TM)  modes  are  possible  only  if  propagation  occurs  parallel  to  the 
principal  axes  of  the  material,  but  that  coupled  TE-  and  TM-like  modes 
occur  for  propagation  at  small  angles  with  respect  to  these  axes. 

In  section  2,  the  formalism  is  given  for  determining  the  modes  for 
propagation  at  an  angle  with  respect  to  one  principal  axis  (y‘)  of  the 
material  and  parallel  to  another  (x'). 

In  section  3,  results  are  obtained  specifically  for  0. 53-pm 
radiation  (doubled  1.06-pm  radiation  from  a neodymium :yttrium  aluminum 
garnet  (NdrYAG)  laser)  in  lithium  niobate  (LiNbOj ) . There,  TE-  and 
TM-like  modes  are  obtained  for  various  thicknesses  of  a slab  backed  on 
one  face  by  a perfect  conductor  and  on  the  other  by  free  space.  For 
some  typical  examples,  the  Poynting  vector,  Sy,  indicating  the  magnitude 
of  power  flow,  is  plotted  versus  distance  into  the  slab  from  the  surface 
to  indicate  the  concentration  of  the  power.  It  is  expected  that  these 
data  when  correlated  with  similar  calculations  for  the  SAW's  will 
determine  more  precisely  the  regions  where  the  interaction  is  strongest. 
These  results  will  help  clarify  the  dependence  of  the  observed  phenomena 
on  polarization  and  the  angle  of  propagation.  Computer  programs  have 
been  written  to  perform  these  calculations  for  the  general  case. 


2.  THEORY 

In  the  cgs  system  of  units,  Maxwell's  equations  in  the  absence  of 
sources  are 
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and  u,  the  permeability,  and  e,  the  dielectric  constant,  may  be  tensors. 
Considering  the  quantities  E,  B,  H,  and  D as  column  vectors,  7X,  y , and 
t as  tensors,  and  V*  as  a row  vector,  we  have,  for  example. 
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Thus,  one  may  easily  see  that  the  product  of  V-  and  VX  gives 
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and  similarly  for  dzx  and  dXy.  Assuming  that  the  functions  (column 
vectors)  that  are  to  be  operated  on  by  V*  (Vx)  are  well  behaved,  the 
partial  derivatives  may  be  interchanged  and 

V • (VX)  = 0 . (7) 

Operating  with  V-  on  equations  (la)  and  (lb),  one  obtains 

V-B  * V-D  = 0 (8) 

as  a natural  consequence  of  Maxwell's  equations  with  no  sources. 

To  solve  Maxwell's  equations  as  a superposition  of  plane  waves  of 
the  form 


exp(ik*r  - i<i)t)  , 


(9) 
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equations  (la)  and  (lb)  become 


kXE  = k pH, 
o 

kXH  = -k  eE, 
o 


where 


k = w/c  , 
o 

0)  is  the  angular  frequency,  and  c is  the  velocity  of  light. 
The  most  general  equation  for  E becomes 

UE  = u, 

where 

U = (kX)vi-1(kX)  + kge  . 

Analogous  to  equation  (4) , 
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Since  equation  (12)  represents  a set  of  three  homogeneous  equations  in 
the  Ex,  Ey,  and  Ez  components  of  E,  solutions  are  obtained  when  the 
determinant  of  coefficients  vanishes. 


det(U)  = 0 . 


(15) 


In  this  report,  we  consider  only  the  case  where  y is  a scalar  and 
equation  (15)  becomes  (multiplying  through  by  y ) 


det  £(kx)  (kX) 

where  from  equation  (14)  one  obtains 
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The  determinant  given  by  equation  (16)  gives  relationships  among  the 
components  kx,  ky,  and  kz  of  k and  w (through  k ) in  the  medium,  and  any 
boundary  conditions  will  further  limit  their  allowed  values.  Once 
values  for  the  components  of  k are  known,  then  all  the  components  of  E 
and  H may  be  solved  in  terms  of  any  one  of  them  by  equations  (la)  and 
(lb)  . 


Taking  the  tensor  components  of  e to  be  .,  where  e.^  = e , then 
equation  (16)  gives  13  "l 

[~(ky  * kz)  * ko  Uc"][-(kz  + kx)  * ko  kc“][-(kK  + ky)  * ko 


[- 

(ky  + kl) 

+ ko  ^nj 
1 

kykz  + ko  P£23 

[- 

fkz  + kx) 

+ k2  ue22J 

kzkx  + ko  pei3 

[-< 

^2  + kA 

x V 

i r 

+ k2  P£ 3 3J 

k k + k2  \1C\0 
x y 0 

- 

*![Vz  + k„2  W2s][> 


k k + k2 


Z X 


’]  [kxky  + ko  ^12]  = 0 


(18) 


This  equation  gives  the  most  general  relationship  between  the  components 
of  k under  the  assumptions  of  an  exponential  space-time  variation  given 
by  equation  (9)  and  of  scalar  permeability  v . 

We  further  limit  the  considerations  in  this  report  to  no  propaga- 
tion in  the  z direction  (kz  = 0),  and  £12  = el3  = 0* 

Equation  (18),  where  k = t\2  = el3  = 0,  reduces  to 
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where 
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(19) 
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The  TE  and  TM  modes  exist  in  the  limit  of  vanishing  £23  (resulting  in  a 
diagonal  dielectric  tensor) . The  TE  mode  for  diagonal  £ is  a solution 
where  Ex=Ey  = 0;  hence,  E = Ez,  which  is  transverse  to  the  direction  of 
propagation  and  experiences  only  the  £33  dielectric  constant  as  shown  in 
equation  (20a) . The  TM  mode  for  diagonal  £ is  a solution  where 
Hx=H  =0;  hence,  H = Hz,  but  E lies  only  in  the  x-y  plane  and  therefore 
experiences  the  combination  of  £n  and  £22  dielectric  constants  shown  in 
equation  (20b) . 

When  £23  t 0,  as  is  considered  below,  these  two  modes  are  coupled  as 
shown  in  equation  (19).  The  specific  problem  that  we  solve  is  a bound 
wave  propagating  parallel  to  the  surface  of  an  infinite  dielectric  slab 
surfaced  with  a perfect  conductor  on  one  face  as  shown  in  figure  1. 


Figure  1.  Representation  of  system  to  be  solved 
for  bound  electromagnetic  radiation  propagating 
in  x-y  plane;  boundaries  only  in  x direction 
are  considered,  with  slab  being  infinite  in 
both  y and  z directions. 


If  the  principal  axes  of  the  dielectric  are  parallel  to  the  x-axis 
but  rotated  about  x through  an  angle  0q  as  shown  in  figure  1,  then  we 
have 
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£11  = 

£1  , 

(21a) 

e22  = 

e2  cos2  eQ 

+ e3  sin2  6o  , 

(21b) 

e33  = 

e3  cos2  e0 

+ e2  sin2  eQ  , 

(21c) 

e23  “ 

(£2  - £3) 

cos  6 0 sin  0 q , 

(21d) 

where  ei  , e2  , and  £3  are  the  principal  dielectric  constants  of  the 
material.  Alternatively,  we  may  consider  - 60  as  the  angle  at  which  the 
electromagnetic  wave  is  launched  with  respect  to  the  e2  principal  axis 
of  the  material. 

In  LiNbOj,  we  choose  z'  as  the  c-axis.  Since  p « 1,  we  have  (ei)*5  = 
( S-2)*3  = 2.34  and  ( e 3 J*3  ■ 2.24  for  0.53-pm  radiation  at  25  C and  (e,)^  = 
( e2 J*5  = 2.24  and  (e3)lj  = 2.15  for  1.06-pm  radiation  at  25  C. 

2.1  Inside  Dielectric  Slab 


equation  (19) , 


Inside  the  slab,  the  solutions  for  k in  terms  of  k are,  from 


2 o 2 2 

kX  = ailC*  kx  =•  k2  r 


(22a) 


where 


2 2 
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k2  - k2 
E M 

1 1 4k23 
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(k2  - k2\2 

V E M/ 

(22b) 


2 2 4 

and  where  kg,  kM,  and  k23  are  given  in  terms  of  k by  equations  (20a,  b, 
c) , respectively.  From  equations  (la)  and  (lb) , one  has 
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(23a) 

(23b) 

(23c) 


10 


Finally,  from  equations  (23d)  and  (24f) , 

s/(e  llko) 


E = -k  H , 
X y z 


(24g) 


Thus,  equations  (24a)  to  (24g)  give  the  six  components  of  E and  H in 
terms  of  the  arbitrary  constants  ulf  u2,  U3,  and  U4,  which  are  to  be 
determined  by  the  boundary  conditions  up  to  a constant  factor. 

2.2  Outside  Dielectric  Slab 

Outside  the  dielectric  slab,  we  have,  in  addition  to  the 
previously  assigned  conditions,  = e22  = e33=  y = 1 and  e23  = 0. 
From  equations  (20a,  b,  c) , which  apply  since  they  are  perfectly 
general,  k^  = 0;  and  the  solutions  of  equation  (19)  split  into  the  TE 
and  TM  modes,  which  give  the  same  solution,  that  is. 


k2  = -k2  + k2 

x y o 


(25) 


We  are  interested  in  a bound  wave,  so  we  choose 


so  that 


k = iT 
x 


r2  = k2  - k2 

y o 


(26) 

(27) 


and  we  look  for  solutions  for  positive  T.  The  choice  of  sign  given  by 
equation  (26)  gives  x variations  of  the  form 


-rx 
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which  die  in  the  positive  x direction  as  shown  in  figure  1 if  T is 
positive. 


The  spatial  solution  outside  the  slab  may  be  written  for  as 


Ez  = -u5e 


-u,e-rxeikyy  . 


(29a) 


Following  the  procedure  for  solving  the  various  equations  starting  from 
equation  (23a) , we  have 


i = r*E  , 
x k z 
o 


H = -i  E 
y k z 

o 


lim  1 /_ 

e23-K)  k2e„,  V 


O 23 


2 2\ 

k_  + k.  -I  r )e2 


2 

‘o  y 


(29b) 


(29c) 


(29d) 


But  since  T2  is  related  to  k£  according  to  equation  (27) , one  gets  an 
indeterminant  form  for  Ey  Thus,  we  may  define  E better  by  choosing  H 
to  have  the  spatial  variation  y 


H = iu  e 
z 6 


r-vv/^.)  , 


(30a) 


and  we  obtain  from  equations  (23d,  e) 


E = - -*H  , 

x k z 
o 


(30b) 


E = i £-  H 
y kQ  z 


(30c) 


These  equations  are  independent  of,  but  analogous  to,  equations  (29a,  b, 
c) . 
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2. 3 Matching  Boundary  Conditions 


We  set  E tangential  and  H tangential  to  be  continuous  across 
the  dielectric-air  boundary  at  x = 0,  and  we  set  E tangential  zero  at  the 
metal  boundary  at  x = -L,  where  L is  the  dielectric  slab  thickness. 


to  get 


From  equations  (24a)  and  (29a),  we  set  E^  continuous  at  x - 0 


U1  + e23U3  = -u5  ? 


(31a) 


from  equations  (24d)  and  (29c),  we  set  continuous  at  x = 0 to  get 


u2  + e23u4  = yFu5  ' 


(31b) 


from  equations  (24e)  and  (30a,  c),  we  set  Ey  continuous  at  x = 0 to  get 


(kl  ‘ Ke)U1  + e23(k2  - Ke)U3 


-E23ru6 


(31c) 


from  equations  (24f)  and  (30a) , we  set  H continuous  at  x - 0 to  get 


(k2  - kjV  + e (k2  - kAu  = - ~ k^  u 

\ 1 E/  2 23V  2 E/  4 e22  M 6 


Old) 


from  equation  (24a) , we  set  = 0 at  x = -L  to  get 


Uj  cos  kjL  - u2 


sin  k^L  sin  k2L 

+ e23u3  cos  k2L  - E23u4  £ — = 0 ' (31e) 


kl 
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For  a solution,  the  following  determinant  of  the  coefficients 
must  vanish: 


k2  - k£ 
2 E 


sin  k^  L 


- 0 kS/E 22  ’ 0 <32) 


sin  k L / .sin  k 

(\2  - kE2) cos  k,L  -(kf  - k£)— irr  (k*1 ' k*)cos  V 


From  equations  (20a,  b,  c) , (22b),  and  (27),  ^each  of  the 
quantities  in  the  determinant  is  expressed  in  terms  of  k . Thus,  the 
solution  of  equation  (32)  gives  the  allowed  values  for  ky. 

In  the  limit  e23  -0  (eQ  = 0 for  propagation  perpendicular  to 
the  c-axis  in  LiNb03),  equations  (31a  to  f)  split  into  two  independent 
sets  of  three  equations.  In  general, 


k2  - k2 
K1  *E 


so  that,  in  the  limit  e23  °» 


2 2 

<kM 


£22(kI  * kE) 


kl  + kE  ' 

k2  •”  kM  ' 
2 2 
kl  -kE  , 


One  solution  is  the  TE  solution  obtained  by  striking  columns  3,  4,  and  6 
in  equation  (32)  and  thereby  getting 


The  general  solutions  for  equation  (32)  are  given  by 


The  second  term  vanishes  as  e23  -*■  0,  since  k|  -*■  k|  , and  the  second  and 

third  factors  in  the  first  term  give,  respectively,  the  TE  and  TM 
solutions  given  by  equations  (36)  and  (38)  (k2  in  this  limit)  . 

2. 4 Transverse  Electric-like  Modes 

We  consider  cases  for  which  0q  is  small  so  that  from 
equation  (21d)  e23  <<  1.  Thus,  we  refer  to  TE-like  modes  as  those 
solutions  of  equation  (40)  where 


tan  k L 


(41) 

(42) 

(43a) 

(43b) 

(43c) 


(43d) 


.ILi 


e22  u5 


In  the  limit  of  small  e23,  equations  (34a,  b,  c)  show  that  u^,  U3,  and 

all  approach  zero,  giving  a pure  TE  mode  for  which  k2  -*■  kg  . 


2. 5 Transverse  Magnetic-like  Modes 

Analogously,  TM-like  modes  are  the  solutions  of  equation  (40) 

for  which 


k2  tan  k,L 


Te22  k2 


1 . 


(44) 


For  these  modes,  we  solve  for  the  u^  in  terms  of  u£  to  get 


^ _ 


■..(‘i  - *2)" 


r - 


tan  k,L 
M 1 


'22 


1 + yr 


tan  kjL 


-i-l 


‘l  J 


(45) 


Then  Uj/Ug  for  1 < i < 4 are  given  by  equations  (43a,  b,  c,  d)  when 
formally  multiplied  by  Ug/Ug.  In  this  case,  in  the  limit  of  vanishing 
e23,  equations  (45)  and  (43a,  b)  show  that  u5 , u3 , and  u2  approach  zero, 
giving  a pure  TM  mode. 

Since  we  consider  cases  for  which  Gjj  > e33  (see  discussion 

2 

following  eq  (21a,  b,  c,  d) ) , there  may  be  solutions  for  which  kM  > °' 


but  k^  < 0 if  ye33J<Q  - ky  * In  these  regions,  k2  > 0 and 

k j < 0,  in  which  cases  kj  = -i|kj|  is  taken  with 


tan  kj  L tanh  ^ | kj  j 


(46) 


Thus,  TM-like  modes  may  appear  for  thickness  L where  TE-like 
modes  are  below  cutoff.  For  pure  imaginary  kJf  no  solutions  for 
equation  (41)  are  possible;  thus,  no  TE-like  modes  are  possible  for 
these  values  of  k2. 

y 
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3.  RESULTS 


Calculations  were  performed  to  determine  the  Poynting  vector,  Sy,  in 
the  direction  of  propagation  and  the  various  field  components  associated 
with  the  lowest  order  modes.  The  ordinary  and  extraordinary  indices  of 
refraction  were  taken  as  nQ  = 2.34  and  n = 2.24,  respectively,  which 
correspond  to  0. 53-vim  radiation  (AQ  = 0.53  x 10-4  cm)  in  LiNb03  at  25  C. 
The  equations  scale,  so  that  the  results  are  identical  for  constant 
Aq/L,  where  AQ  is  the  free  space  wavelength  of  the  radiation.  Several 
angles  (0O)  were  chosen  for  the  orientation  of  the  dielectric  axes  with 
respect  to  the  direction  of  propagation,  which  was  kept  fixed  along  the 
y-axis  as  shown  in  figure  1.  The  angle  0Q  affects  the  parameters  e^j 
accordir  to  equations  (21a  to  d) , for  which 


e3  = n2  = 5.0176. 
e 


ei  = Co  = n2  = 5.4756  and 
1 z o 


For  small  angles  of  propagation  of  about  0O  = 5 deg  or  less,  the 
results  are  essentially  the  same  as  for  0Q  = 0 deg  in  which  case  the 
axes  are  collinear.  For  0O  = 0 deg*  the  fundamental  TE  and  TM  modes  for 
AQ/L  = 0.53  are  shown  in  figures  2 and  3.  Although  the  total  power 
flows  are  not  normalized,  one  can  distinguish  that 


a.  The  TE  mode  (fig.  2)  has  only  one  component  of  the  E-field  that 
is  along  the  dielectric  z-axis  parallel  to  the  surface  of  the  slab. 


sees"  predominantly  ne,  except  for  the  energy  in 


b.  This  mode 
free  space;  hence,  the  effective  index  of  refraction, 
less  than  n 


ngff  = 2.22553. 


nef f , is  slightly 


c.  The  TM  mode  (fig.  3)  has  two  components  of  E,  one  perpendicular 
to  the  surface  of  the  slab  and  one  parallel  to  the  surface  along  the 
direction  of  propagation  (y) . 


d.  This  mode  "sees"  predominantly  nQ,  except  for  the  energy 
free  space;  hence,  neff  is  slightly  less  than  nQ: 


in 


n , , = 2.3363. 
eff 


e.  The  power  flow  of  the  TE  mode  is  closer  to  the  surface  of  the 
»slab,  and  a greater  percentage  of  its  energy  is  contained  in  free  space. 

In  contrast,  the  power  flow  of  the  TM  mode  is  closer  to  the  conductor. 

f . The  magnitude  of  Ey  in  the  TM  mo^e  relative  to  Ex  gets  larger 
as  the  surface  is  approached  from  within  tl  dielectric.  (This  property 
may  be  significant  in  understanding  the  interaction  of  these  waves  with 
SAN's.) 
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A 


The  modes  depicted  in  figures  2 and  3 can  be  labeled  the  TE0  and  TMq 
modes  in  which,  in  general,  TEm  and  TMm  modes  have  m changes  in  sign  in 
theirEzand  Hz  components,  respectively.  As  L increases,  the  relatively 
little  change  in  the  form  of  the  modes  indicates  that,  within  a fixed 
distance  from  the  surface  of  the  dielectric,  there  is  proportionately 
less  fractional  power  of  the  total  flowing  in  the  dielectric.  For  the 
higher  order  modes,  m large,  a larger  amount  of  power  is  in  the  free 
space  near  the  surface.  As  m grows  larger,  nfiff  1,  k^  -*■  k£,  and 
r ■+  0;  thus,  the  damping  of  the  evanescent  portion  of  the  wave  in  free 
space  is  less  severe. 


(E//TO  SURFACE) 
(n^f  • 2 22553) 


Figure  2.  Electromagnetic  field  components  for 
fundamental  transverse  electric  (TEQ)  mode 
in  slab  of  thickness  L = Xo/0.53  for  0g  * 

0 deg  angle  of  propagation. 
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The  variation  of  rigff  of  the  TE^  mode  for  different  L as  a function 
of  0O  is  shown  in  figure  4.  The  thickness  Lm  corresponds  to  an  "m"  ym 
thickness  slab  for  XQ  = 0.53  vim  (green)  radiation;  thus,  the  infinite 
medium  case  is  rapidly  approached  by  small  thicknesses.  For  practical 
thicknesses  (-1/8  in. — 0.3  cm),  laser  radiation  incident  as  a sheet  beam 
over  the  surface  of  the  dielectric  excites  a very  large  number  of  higher 
order  modes.  The  variation  of  the  field  components  near  the  surface, 
however,  is  very  much  like  the  variations  near  the  surface  depicted  in 
figures  2 and  3. 


POYNTING  VECTOR 


ANISOTROPIC  MEDIUM 


lE^TO  SURFACE) 
(ntf|  • 2 3363) 


Figure  4.  Variation  of  effective  index  of  refraction 
(neff)  of  fundamental  transverse  electric  (TE)  mode 
for  various  slab  thicknesses  (1^,,)  as  function  of 
angle  of  propagation  (9p) ; neff  °f  fundamental  trans- 
verse magnetic  mode  is  almost  independent  of  9q  and 
varies  as  2.3363,  2.3391,  2.3396,  and  2.340  for 
L = L, , L,,  L,,  and  L , respectively. 

1 Z O oo 


Figure  5 shows  a typical  higher  order  mode  propagating  at  an  angle 
of  = 19  de<?*  This  is  a TMj-like  mode  (two  zero  crossings  for  Hz), 
which  shows  the  relative  amounts  of  Hx,  Hy,  and  Ez  components  generated 
at  this  angle.  These  components  are  absent  in  TM  modes  propagating  at 
the  9 q — 0 deg  angle.  The  discontinuity  in  Ex  at  the  surface  is 
clearer  in  this  higher  order  mode  than  in  the  TMQ  mode  in  figure  3. 
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0 0 - 19  DEG 


ANISOTROPIC  MEDIUM 


(n,„  » 2.245751 


Figure  5.  Electromagnetic  field  components  for  transverse 
magnetic  (TM2)-like  mode  in  slab  of  thickness  L = XQ/0.53 
for  = 19  deg  angle  of  propagation. 


4.  CONCLUSIONS 

These  calculations  show  that  many  TE  and  TM  modes  are  allowed  in 
slabs  of  practical  thicknesses.  The  problem  of  forming  linear 
combinations  of  these  modes  to  approximate  a sheet  laser  beam 
propagating  near  the  surfaces  of  thick  slabs  is,  accordingly, 
formidable.  What  is  necessary  for  a qualitative  understanding  of  the 
interaction  between  such  sheet  beams  and  SAW's  may  be  clear,  however, 
from  the  examples  given  here.  If  further  work  is  necessary  for  a more 
quantitative  picture,  the  formalism  derived  here  and  the  computer 
programs  developed  in  this  connection  are  available  from  the  author. 
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